m 



PROCEEDINGS 

°F SCIENCE 



Recombination in nuclear collisions 



o 

(N 

Oh' 
(D 
00 

(N 



Rudy Hwa* 

Institute of Theoretical Science, University of Oregon, Eugene, OR 97403, USA 



E-mail: hwa@uoregon.edu 



> 

o 
p 

o 
o 



X 



Recombination is a hadronization process that converts partons to hadrons at late time, but the 
description has no quantitative significance without some meaningful input on the parton dis- 
tributions at earlier time. Thus observations of particle spectra and correlations have definitive 
implications on the partonic processes at all transverse momenta. After presenting a general re- 
view of the subject at the Workshop, I selected two topics in nuclear collisions for more detailed 
discussion, which are summarized here. One is on the azimuthal anisotropy at low pr due to hard 
or semihard scattering of partons that create ridges with or without triggers. The ridge particles 
are the products of recombination of thermal partons enhanced by the energy loss of hard or semi- 
hard partons. Their (j) dependence at midrapidity can be determined essentially from geometry. 
The other topic is on the scaling behavior in (j) and centrality at high pr, where the hadronization 
process is dominated by TS and SS recombination. The relevant RHIC data are found to have 
the same scaling behavior. But at LHC such scaling is badly broken at pj ^ 10 GeV/c if two-jet 
recombination is important. At the end of this article some comments are made to relate our study 
of the effects of semihard partons to the observation of minijets in the analyses of STAR data. 
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1. Introduction 



In this workshop various approaches to the study of nuclear collisions have been presented 
and compared. I have been asked to present the approach that emphasizes parton recombination. 
Although hadronization by recombination is a process that occurs at late times when the medium 
density is substantially reduced from the early time, the problem of determining what partons 
recombine opens up the study to include the physics at all times. Indeed, while the media created in 
leptonic, hadronic and nuclear collisions may differ drastically initially, the hadronization processes 
in the final stages of their evolutions are essentially universal. Of course, fragmentation is uniquely 
suitable for leptonic processes, but that is because of the absence of neighboring medium partons. 
In a more general approach to hadronization, fragmentation is a subclass of recombination. Thus 
the issue returns to the question of what recombines, and the partonic system that is formed by a 
collision is the real object of investigation. 

At the suggestion of the organizers my talk included topics on the formulation of recombina- 
tion in the beginning many years ago. The valon model was introduced as a link between deep 
inelastic scattering at high and hadronic wave function at low Q^, which in turn made possible 
a description of the recombination function in the time-reversed process. Application was then 
made to the low-pj- process of p + p ^ n + X. Shower parton distributions were derived from the 
fragmentation functions (FF) treated as the result of recombination of shower partons. Such an ap- 
proach led to the possibility of relating the FFs of partons into mesons (D*^) and baryons {D^), thus 
closing a conceptual gap in the relationship of the FFs among various species of hadrons produced. 

In nuclear collisions one usually identifies three domains of pT where different approaches are 
separately most suitable: low region {pT < 2 GeV/c) where hydrodynamics has been successful, 
high region {pj > 6 GeV/c) where lower-order calculations in pQCD is valid, and intermediate 
region {2 < pj <6 GeV/c) where recombination dominates over fragmentation. Since the ratio of 
/D^ is small, it has been hard for the fragmentation paradigm to reproduce the large p/n ratio 
at intermediate px in nuclear collisions. Even in pA collisions the Cronin effect is larger for pro- 
duced proton than for produced pion and therefore puts considerable difficulty in the conventional 
interpretation that the effect is due only to initial-state transverse broadening followed by fragmen- 
tation. Recombination can reproduce the data quite readily. Space does not permit me to review 
here all those topics. A summary of them with extensive references to the original articles can be 
found in ijl]]. In the text below I discuss the more recent development in nuclear collisions. I chose 
two topics of particular relevance to this Workshop because they depart from the conventional wis- 
dom. One is on azimuthal dependence at low pj where the usual explanation is that the variation 
of pressure gradient on <p leads to elliptic flow. We consider the effects of semihard scatterings that 
precede the completion of thermalization? The other topic is two-jet recombination at LHC that 
can overwhelm fragmentation for hadrons formed at pj around 10 GeV/c. 

At the end I offer some thoughts on possible common ground with the 2-component model 
of the UW-UTA alliance in their interpretation of the inclusive distribution and auto-correlation 
determined in their analyses of the data from heavy-ion collisions at RHIC [^. 
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2. Azimuthal Dependence at Low pj 

It is conventionally thought that for pj < 2 GeV/c only soft physics is involved and hydrody- 
namics provides a satisfactory description of the elliptical flow coefficient V2{pt) ^ 0- However, 
the theory assumes rapid thermalization, achieved at Tq = 0.6 fm/c. Whether that can be proven is 
a separate issue. What is relevant to the following discussion is the question on whether semihard 
scattering with parton kj less than 3 GeV/c can have any effect that supersedes the hydro result, 
since there are many such scatterings and they take place in T^O.l fm/c. The point is that whereas 
hard scatterings are rare and can be ignored in the treatment of the bulk, there may come a point 
when the transverse momenta kj are low enough such that semihard partons are pervasive and not 
negligible, yet not accounted for by the hydro treatment of a macroscopic medium in local thermal 
equilibrium. A possible empirical basis upon which such questions on the "conventional wisdom" 
can be raised is the implication of the existence of minijets by the two-component analysis of the 
STAR data [§, ||]. Without necessarily endorsing the specifics of the two-component model, we 
examine the possibility that semihard partons emitted from the surface can lead to the observed 
azimuthal anisotropy, even if thermlization cannot be completed until after 1 fm/c. 

Our focus in this section is on ridge production at midrapidity, given the experimental fact that 
the extraction of the (p dependence of associated particles relies on a section of Arj > that char- 
acterizes the ridge (elongation in Arj) [^. Ridge formation at large Arj is a separate problem and 
is to be addressed separately Here our emphasis is in the study of (p dependence in the single- 
particle distribution pi [pj , , ^) at Tj 0, arising from semihard partons without any triggers. Our 
point is that if semihard partons are copiously produced near the surface with a <p anisotropy that 
is calculable from geometry, then they can have nontrivial consequences on p\{pT,<p,b) in a way 
similar to how ridge structure can develop in the particle distribution associated with a high-/7r 
trigger [||,^ [I^. Although the use of ridge in the nomenclature referring to such structure was ini- 
tiated by Putschke ||^, the discovery of the same behavior (narrow in but broad in Ax]) occurred 
earlier in autocorrelation without triggers [pl|]. 

When a semihard scattering occurs deep in the interior of the nuclear overlap, the scattered 
partons get absorbed by the medium and become a part of the bulk. If it occurs near the surface, 
one of them can get out after losing some energy, which in turn can enhance the thermal energy 
of the partons near its trajectory. Recombination of the enhanced thermal partons gives rise to the 
ridge particles locally. Since there can be many semihard scatterings in each event of central or 
mid-central collision, the sum of all such ridge particles produced gives rise to a dependence in 
Pi{PTA,b). Our assumption is that the base upon which the ridge sits is independent. Thus we 
write 

Pi{pT,^,b)=B{pT,b)+R{pT,^,b). (2.1) 

The first term on the right side is not the hydro bulk, which has dependence and is all of 
Pi{pT,<p,b) at pt < 2 GeV/c. We call it Base, which is independent, on the grounds that the 
anisotropic geometry of the overlap at ft > leads to calculable asymmetry of the semihard par- 
tons near the surface at early time and produces the only (p dependence of p\{pT,^,b) at late time 
through the ridge term R{pT ,<p,b).ln that picture it is necessary to demonstrate that the ridge struc- 
ture in the single-particle distribution can be related to the ridge yield in two-particle correlation 
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with a detected high-/7r trigger, since the same physics is involved. In that way we have a unified 
picture without patching together different pieces of physics to distinguish different kinds of ridges. 

For notational simpUcity, let us use p^(0i,02) to denote the ridge part of the two-particle 
distribution with a hard or semihard parton at 0i and several thermal partons at ^2 that can form 



the detected hadron (pion or proton) with the same 02- The last term in (2.1), abbreviated as R{<^), 



gives the ridge distribution whether or not the parton at 0i leads to a trigger. Thus we have 

R{h)^ j d^ipa^uh)- (2.2) 

On the other hand, if events are selected by trigger particles at (^1 , the ridge yield is integrated over 
all associated particles at 02 > so we should have for the yield 

F«(0i)°c|j02pf(0i,02). (2.3) 

Our aim is to give a theoretical description that quantifies the above relations as well as to show 
that they are consistent with the data at hand |[l0|]. 

The basic premise of our approach is that the elliptical nuclear overlap of the initial config- 
uration in the transverse plane at j is sufficient to describe the dependence of the ridge 
distribution without carrying out a detailed hydrodynamical calculation. Our assumption is that 
on average the final hadrons are directed normal to the surface of the initial ellipse. Thus for any 
point on that surface with its normal vector denoted by one can ask what semihard partons at 
01 can contribute to the formation of a ridge particle at This question was addressed earlier in 



an article on the correlated emission model (CEM) [12], where a Gaussian width of a f« 0.33 was 
found to give a good fit of the 0i dependence of the ridge yield [13|. That is, A0 = 02 — 0i cannot 
be too large if there is to be coherence of the enhanced thermal partons to form a detectable ridge 
particle. If A0 is much larger than a, then the effect of the energy loss of the semihard parton is 
spread out over a large angular range due to the radial flow of the thermal partons, too diffused to 
result in effective coalescence. The physics can be simply summarized by a line segment 5(02)^) 
on the surface around 02 through which semihard partons can contribute . For every b, S{<^,b) can 
be analytically determined in terms of elliptic integrals of the second kind [|^], the details of which 
need not be reproduced here. Its dependences on for = and 7r/2 are shown in Fig. 1. For 
mid-central collision (ft 1 in units of Ra) where the initial elUpse is neither too narrow nor too 
broad, there is a large difference between emission at « and that at « 7r/2. That sounds like 
the result of hydro calculation, but no hydrodynamics has been applied. 

Since the probability of producing the ridge particles depends on the density of semihard 
partons that give rise to the enhanced thermal partons, we thus have 

R{pTA.b)ocS{<^,b). (2.4) 

The proportionality factor involves the nuclear density and an exponential dependence on pj. Since 
S{^,b) prescribes the only dependence of pi (/jj- , , ft) , we can determine the strength of the ridge, 
relative to the base, by fitting V2{pT,b) jlOl]. That is shown in Fig. 2. The ridge yield, on account 



of (|2.3D, also depends on the surface segment according as 

F'^(0i)oc5(0i,ft). (2.5) 
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Figure 1: Surface segment S{(p,b) vs normalized impact parameter b in units of for ^ = (solid line) 
and (j) — 7t/2 (dashed line). 
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Figure 2: Elliptic flow coefficient V2{pT,b) vs b for pj — l.2\ GeV/c (solid) and 1.71 GeV/c (dashed). Data 
are from 0. 



Fig. 3 shows how well the data are reproduced by ( |2.5[ ) with normalization adjusted to render a 
good overall fit of the 0-5% data [Fig. 2(a)]. Both the normalization of the 20-60% data [Fig. 2(b)] 
and the shapes for both centralities are well reproduced by the calculated results without further 
adjustment. This demonstrates that the triggered and untriggered ridges are related by the same 
surface segment S{<p,b), without requiring rapid thermalization or pressure gradients. 



3. Scaling Behavior at High pr 

At high pt it is necessary to take into account hard partons traversing the nuclear medium. 
PHENIX has data on Raa that are separated into bins of pt,({> and Np^n [Q- The dependencies 



5 



Recombination in nuclear collisions 



Rudy Hwa 



0.2 



0.15 



DC 

> 



0.1 



0.05 



Q 



1 1 1 1 

(a) • 0-5% 


1 1 1 1 

(b) ■ 20-60% 


STAR 




1 1 1 1 





20 40 60 80 20 40 60 80 



Figure 3: Ridge yield per trigger vs trigger angle 0s relative to reaction plane for (a) 0-5% and (b) 20-60% 
centralities. Solid lines are calculated results given in [ |lo[ ] based on (2.5), and the data are from |jl3|]. 



on the three variables are so complicated that there is a need for organizational simplification. 
Obviously, ^ and A'part (or b, as we prefer to use) are related by geometry. We find that an average 
dynamical path length can be defined, in terms of which Raa has universal scaling behavior that 



can in turn be reproduced in our theoretical calculation [15]. 

Consider a collision at b with a high-/7r particle detected at 0. Suppose that a hard parton is 
produced at the position XQ,yQ in the initial elliptical overlap with a transverse momentum k directed 
at . We consider only the transverse plane at mid-rapidity. Let us define a geometrical path length 
by 



^{xQ,yQ,(^,b) 



ti(xo,yo,<l),h) 



dtD[x{t),y{t)], 



(3.1) 



where t is just a label to mark the trajectory from the point of origination to the exit point ?i. D is 
the density of the medium that the hard parton traverses, and can be calculated from the almond- 
shaped overlap region. Clearly, if the nuclear medium has hypothetically low density, it is sensible 
to regard i as zero even though ti ^ 0. We do not consider medium expansion for two reasons: (a) 
if we did, it would be necessary to involve hydrodynamics that is an additional theoretical input 
we want to avoid, and (b) an expansion would lower D but lengthens t with compensating effects. 
There are subtleties in mapping the almond to the elliptical boundaries, especially when (xojjo) 



is not too far from the surface [15]. For long transit time, which is less important, longitudinal 
expansion should be taken into consideration, and we rely on the dynamical scaling behavior found 



in [16] to regard our procedure as being insensitive to that expansion. Basically (3.1 ) is precise and 
essentially geometrical. 

To take the dynamics of momentum degradation into account, we define the dynamical path 
length 



§ = 7£(xo,yo,0,^) 



(3.2) 
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Figure 4: Experimental data on Raa for pT=^-5 (left) and 7-8 GeV/c (right) from Ref. [|l4[ are plotted 
against the scaling variable ^ . The values are defined hy ^ ~ nn/2A with n being 1,3,.. from left to right, 
all having the same symbol for each centrality. The solid lines are the same as the solid lines in Fig. 5 that 
represent the results from theoretical calculation. 

with 7 to be determined, and the average is 

ft) = y j dxQdyQe{xo,yo,(l>,b)Q{xo,yo,b), (3.3) 

where Q{xo,yo,b) is the probability of producing a hard parton at xo,yo and can be determined 
from the thickness functions of the colliding nuclei at b. The significance of considering f (0,^) 
is that for every pair of {<p,b), one can calculate | and therefore /?aa(|)> given the data point at 
{(j),b). What we have found is that for every bin of pr the data on Raa exhibit universal behavior in 
^, independent of the variations of {<p,b) for each ^. That is shown in Fig. 4, which is a property 
of the data independent of theoretical input. Only two bins of pr are shown, but they exemplify 
the dependence on f for 10 bins in the interval 2 < pj < 10 GeV/c. A specific value of 7 has 
been used in the plots, but that should not affect the existence of the scaling behavior, since 7 only 
determines the horizontal scale. Obviously, the scaling behavior provides considerable economy 
in the presentation of the data, and makes possible an efficient comparison with the results of 
theoretical calculation. 

For pt > 2 GeV/c it is necessary to consider the shower partons, which in turn depend on 
the hard parton distribution Fi{q,^) at the surface with momentum q. It is related to the parton 
momentum k at the point of creation by 

Fi{q,^) = I dkkfi{k)G{k,q,^), (3.4) 

where fj{k) is the k distribution of parton type /, and G{k,q,^) is the degradation function 

G{k,q,^)=q5{q-ke-^). (3.5) 

The <p dependence is completely imbedded in <^ . In the usual recombination formalism the shower 
parton distribution can be determined from Fi{q,^). For the thermal partons the surface segment 
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Figure 5: Theoretical result showing scaling behavior of Raa at pr = 4.5 GeV/c (left) and pT = 7.5 GeV/c 
(right) for six values of centrality. The points for different values between and n/2 have the same symbol 
when the centralities are the same, with increasing from left to right with n= 1,3, •••,11. The solid lines 
are the best fits; the dashed line is an exponential fit by (3.7). 

function S{^,b) in Fig. 1 is needed to incorporate the dependence. With those ingredients the TS 
and SS recombination components can be calculated for the single -particle distribution pi{pT,^) 
and then for Raa{^) [p3|]. The results are shown as points in Fig. 5 for six values of centrality 
(denoted by c) and six values of (p (defined by nn/24). All 36 points lie on a universal curve 
for each bin of pj, when plotted against i§. The value of y is adjusted to fit Raa{pt, ^,Npan) for 
4 < pt < 5 GeV/c at various values of (p and A'pan; it is found to be 

7 = 0.11. (3.6) 

Thus one set of data has been used to fix the degree of momentum degradation in the nuclear 
medium; all other results are obtained without further adjustment. 

The solid lines in Fig. 5 are analytic curves to represent the calculated points. Those curves are 
transported to Fig. 4 to facilitate the comparison between theory and experiment. The agreement 
is evidently very good. Thus the recombination model satisfactorily reproduces all the data points 
for all Pt > 2 GeV/c and all <p and c with one adjustable parameter 7, whose implication on 
momentum degradation can be more simply described by the dashed line in the left panel of Fig. 
5. It is a straight-line fit of /?AA(f )> giving 

^5A(^")|p,=4.5 = ^'^P(-2-6f), (3.7) 

describing an exponential suppression at increasing f , which is physically reasonable, and is de- 
pendent only on the average path length, independent of separate values of (p and b, which is not 
surprising, yet nevertheless a new finding. It is hard to relate the numerical value of 7 to other 
physical quantities that we are familiar with. 

In extending our consideration to LHC we find that with all the physical mechanisms remain- 
ing the same, the scaling feature in ^ is preserved, as shown in Fig. 6(a) in a linear plot for pj = 10 
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Figure 6: (Left) Scaling behavior of Raa at LHC for pT=^Q GeV/c and negligible 2-jet recombination. The 
symbols are the same as in Fig. 5 for various values of c and 0. (Right) Non-scaling behavior of Raa at LHC 
for pj=\Q GeV/c when 2-jet recombination is dominant. 



GeV/c. However, the hadronization process need not be the same. At LHC the density of hard 
partons created is so high that the shower partons from neighboring jets can recombine to form 
hadrons. We do not mean large xt, for which parton density is low whether at RHIC or LHC. 
At fixed pt the parton density increases with y^, so it is possible that at some energy the shower 
cones of adjacent jets may overlap. In that case the recombination of shower partons arising from 
parallel jets can produce hadrons copiously, since the hard-parton momenta, and k2, can be less 
than pj, as opposed to being greater than pj for fragmentation. In Fig. 6(b) we show the result for 
Pt = 10 GeV/c at LHC for a particular choice of the overlap parameter T in two-jet recombination 
[]T5|]. There are two striking features in that figure. One is that the scaling behavior in Fig. 6(a) is 
badly broken. The other is that the magnitude of Raa is extraordinarily large, with the possibility 
of exceeding 1. The observation of either feature would be a clear signature of significant depar- 
ture from RHIC physics, the latter would stand out distinctly from other predictions. The origin of 
scaling violation and large Raa is in the definiton of Raa, where the normalization factor involves 
NcoW for single jets. When 2-jet recombination becomes important, normalizing Raa by N^oW ^i^^ 
restore f scaling and lower Raa- We add that Raa for p production at LHC can be even larger due 



to the recombination of three shower partons [|17[]. 

It should be pointed out that if the features in Fig. 6(b) are not found at LHC, it does not imply 
the failure of the recombination model. It just means that the overlap of neighboring jet cones is 
not of high probability to enhance 2-jet recombination. The observation of any features between 
Fig. 6(a) and (b) would provide some empirical ground for quantifying the recombination function 
of adjacent jets. 



4. Comparative Comments 

Since the aim of this Workshop is to exchange ideas on different approaches to RHIC physics, 
it seems appropriate to end this talk with some comments on the similarities and differences among 
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some of the models. In particular, there seems to exist some common ground between the point 
of view advanced in Sec. 2 and the two-component model advocated by the UW-UTA alliance 
That model consists of a soft S and a hard H component with the most notable observation that 
in Au-Au collisions H has a strong enhancement at small transverse rapidity yt. It is interpreted 
as being due to minijets which peak at around = 2.7, corresponding to pj ^ \ GeV/c. Putting 
aside any reservation about the procedure of separating the two components, the observation of a 
pronounced peak in the yt-jt correlation and the necessity of including a 2D Gaussian peak in order 
to fit the autocorrelation data on t/a-^a in heavy-ion collisions provide some empirical evidence 
for the existence of minijets [|2|]. It is claimed that the hadron spectra of the hard component reveal 
the dominance of parton fragmentation all the way down to below ~ 1 GeV/c, and that the 
abundance of surviving minijets conflicts with near-ideal hydrodynamics [|5|]. 

In Sec. 2 we have considered the possible effect of semihard partons on anisotropy before 
thermal equilibrium is established. Since semihard scattering is not rare for kj <?> GeV/c, there 
is a possibility that there may be a correspondence between those semihard partons and minijets. 
The hadronization mechanisms may be quite different, ours being the recombination of enhanced 
thermal partons, while the 2-component model stresses fragmentation. Our ridge component in 
the single-particle distribution contains the dependence, corresponding to the hard component 
that has the low-pr enhancement. AutocoiTclation exhibits 0a dependence after integation over 

It would be of interest to know the dependence on ^ before integration (admittedly a bit hard 
to achieve), since that would give a possible hint of the dependence of S{<p,b), the boundaries 
of which are shown in Fig. 1, but it is totally known analytically [p^. The ridge yield has the 
curious property that, when the trigger is along the direction of the reaction plane {(pg = 0), the 
yield increases as the collision becomes less central. That is a property shown by the STAR data in 
Fig. 3 and corresponds to the solid line in Fig. 1. Since S{<p,b) is the surface segment for semihard 



partons whether they lead to trigger [cf. ( |2.5[ )] or not [cf. ([2.4|)], it would be good to verify its 
relevance to both the triggered events [ pj| ] and the untriggered autocorrelation. 

There are numerous questions in my mind about the 2-component model when the collision 
energy is raised to those at LHC. As the semihard parton density becomes high, the recombination 
of adjacent jets will at some point become unavoidable. The separation of soft and hard com- 
ponents will also become more difficult, especially if the hard component is so identified with 
fragmentation, as done at RHIC [^. But that will be a problem for all theoretical models that 
exclude recombination. Of course, recombination is only one of many paradigms that has not yet 
been proved to be universal. I encourage the organizers to reconvene this Workshop in two or three 
years and review whether the RHIC paradigms will have significantly changed at LHC. 
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